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ABSTRACT 

Context. High-resolution, high-sensitivity multifrequency radio images of supernova remnants (SNRs) are essential in advancing 
the understanding of both the global SNR dynamics and particle acceleration mechanisms. 

Aims. In this paper we report on a new study of the SNR Puppis A based on VbA observations at 1425 MHz; the improvement 
represents a factor of two in angular resolution and almost ten times in sensitivity compared to the best previous image of 
Puppis A. This new image is used to compare with re-processed 327 MHz data and ROSAT and Chandra images to investigate 
morphological and spectral characteristics. 

Methods. The observations were carried out with the VLA in the DnC and CnB configurations in 2004. After combining with 
single-dish data from Parkes, an angular resolution of 34" x 16" and an rms noise of 0.5 mjy beam -1 , were achieved. Archival 
VLA data at 327 MHz were also processed. The spectral index distribution was then determined by a direct comparison of 
the homogenized data at both 327 and 1425 MHz. In addition, to identify different spectral components, tomographic spectral 
analysis was performed. 

Results. The new 1425 MHz radio image reveals a highly structured border encircling a diffuse, featureless interior. In particular, 
the northern half of Puppis A displays a complex structure along the periphery, consisting of short arcs resembling "wave-like" 
features. These are oriented essentially perpendicular to the shock front on the NE side, but are tangential to the shock on the 
NW side. A remarkable correspondence between such "wave-like" features and spectral changes is observed. On the other hand, 
the brightest radio features (located to the E of the SNR and also detected in X-rays) have no counterpart in the spectral index 
distribution. Based on a uniform compilation of integrated flux densities between 19 and 8400 MHz, a global spectral index 
a — —0.52 ± 0.03 (S oc v a ) has been determined. The new 1425 MHz image of Puppis A was compared with the ROSAT X-ray 
image of the entire SNR and with the resolved arcsec Chandra image of the bright E region. There is good overall agreement 
between the radio and high resolution X-ray images. At the arcsec level, however, the agreement of the radio and X-ray images 
is less striking. A noticeable similarity is found between Puppis A, with its E and W extensions, and analogous morphological 
features observed in the SNR W50. This suggests that Puppis A could be another case of a SNR shaped by the action of 
energetic jets. 
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1. Introduction 

Sensitive, high-resolution studies of the radio emitting 
plasma associated with supernova remnants (SNRs) have 
proven to be an excellent tool for recognizing morpho- 
logical characteristics, delineating the location of contact 
discontinuities in the fluid, and spatially resolving fila- 
mentary emission that may indicate expansion into non- 
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uniform environments in the interstellar medium (ISM). 
These high quality radio images are also helpful in ad- 
dressing questions concerning particle acceleration mech- 
anisms based on determinations of the spatial variations 
in the radio spectral index. Moreover, the detailed combi- 
nation of radio and X-ray data not only complements the 
analysis for particle acceleration, but is also a valuable 
tool for studying the interaction of the blast wave with 
inhomogencitics in the surrounding ISM. 

Dubner et al. (1991) have presented an earlier VLA 
image of the SNR Puppis A at 1515 MHz based on ob- 
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servations from a combination of 26 different pointings 
observed with the VLA 1 in the DnC array. The angular 
resolution of this image is 77" x 43" with a noise of ~ 3 
mjy beam" 1 . This image shows some agreement between 
the bright radio features and the X-ray emission detected 
by Einstein X-ray Observatory. However, the angular res- 
olution and sensitivity of this previous VLA image is not 
high enough to allow a detailed comparison with improved 
X-ray images, e.g. that of Petre et al. (1996) obtained us- 
ing ROSAT, or the detailed high resolution Chandra study 
of the eastern region carried out by Hwang et al. (2005). 
Because of advances in both instrumentation and imaging 
algorithms for the production of sensitive radio mosaics of 
large regions, a new set of 20 cm radio observations have 
been carried out using the VLA. 

In this paper, we present a new 1425 MHz image of 
Puppis A generated from VLA DnC and CnB configura- 
tion data sets obtained in 2003 and 2004. The improved 
sensitivity and resolution have led to the discovery of new 
morphological features near the outer shock front and at 
the E side of the SNR shell. From the combination of the 
radio observations at 1425 MHz with re-processed data 
at 327 MHz from the VLA archives, we have investigated 
the spatial spectral index distribution across the remnant. 
Based on the 1425 MHz image, we have also performed a 
detailed study of the correlation between radio and X-ray 
emission, focusing special attention on the complex region 
to the E. 

1.1. General multiwavelength properties of the SNR 
Puppis A 

Puppis A is the remnant of a galactic supernova for which 
it is possible to observe ejecta from the original explo- 
sion that occurred ~ 3700 yr ago (Winkler et al. 1988). 
In addition, signatures of the interaction with the circum- 
stellar/interstcllar material can be observed. Based on an 
analysis of both X-ray and optical spectra, Puppis A can 
be classified as belonging to the small class of galactic 
oxygen-rich SNRs (Canizarcs & Winkler 1981; Winkler & 
Kirshncr 1985) in which fast moving filaments evidence 
nuclear processes that occurred in the progenitor star. 

The radio emission of Puppis A is of synchrotron origin 
and takes the form of a highly distorted clumpy shell, with 
the brightest section along the E border. Optical emission 
is present as a collection of compact condensations rather 
than delicate filaments, and shows little relation to the 
radio features (Milne et al. 1983). 

Previous multiwavelength studies of the surroundings 
of Puppis A suggested the existence of dense interstellar 
clouds along the E-NE border of the SNR (Petre et al. 
1982; Teske & Petre 1987; Arendt et al. 1990, 1991). This 
was later supported by atomic and molecular studies of 
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the surrounding interstellar medium; HI and 12 CO J: 1-0 
line observations have revealed the existence of a chain of 
interstellar clouds concentrated along the NE periphery of 
the remnant (Dubner & Arnal 1988; Reynoso et al. 1995). 
From these studies, a kinematical distance of 2.2 ± 0.3 
kpc has been derived for Puppis A. 

The X-ray remnant consists of non-uniform emission 
that includes both extended features and knots (Petre 
et al. 1982), the most conspicuous of which is the bright 
knot located at the E limb (called by those authors "the 
bright eastern knot", or "BEK"). This feature coincides 
with an "indentation" in the radio morphology, suggest- 
ing an encounter between the SN shock and an external 
inhomogencity. Based on an OH 1667 MHz line study, 
Wocrmann et al. (2000) suggested that the shock front 
expanding to the E has completely engulfed a cloud giv- 
ing rise to the X-ray knot, leaving no observable unshocked 
remainder of the cloud. Recently, Hwang et al. (2005) pre- 
sented Chandra X-ray images and spectral analysis of the 
region around the BEK. These images reveal details of 
substructure in this part of Puppis A. The spectral analy- 
sis shows a temperature gradient rising behind the shock 
front, suggesting possible heating by reflected shocks due 
to the interaction with denser gas. 

An unresolved X-ray source has been detected near the 
centre of Puppis A (RX J0822-4300, Pavlov et al. 1999, 
and references therein) . This source may well be the com- 
pact stellar remnant formed in the SN event. However, no 
associated pulsar wind nebula has been detected (Gaenslcr 
et al. 2000). 

2. Observations 

2.1. Re-processing of radio continuum data at 327 
MHz 

327 MHz data acquired with the VLA in the CnB configu- 
ration in 1988 (Dubner et al. 1991) were re-processed. The 
data were flagged and calibrated using the AIPS software 
package. Further processing was done in AIPS++. The W- 
Projection imaging algorithm (Cornwcll et al. 2003) was 
used to correct for the non-coplanarity of the VLA array 
at this frequency (Cornwell & Pcrlcy 1992). We made use 
of multiscale clean to deconvolve the image. The decon- 
volution was done in iterations, with the visibility data 
being self-calibrated in phase after each iteration. 

To recover the very low spatial frequencies that the 
VLA observations miss, we made use of single-dish data 
extracted from the 408 MHz All-Sky Continuum Survey 
(Haslam et al. 1982). We scaled the flux density of the sin- 
gle dish image with a spectral index a = —0.6, which we 
derived from comparing the flux densities at 86 and 408 
MHz (Milne et al. 1993). The resultant single-dish image 
was combined with our interferometric image. The com- 
bination is based on a method referred to as "feathering" . 
This method involves Fourier transformation of both the 
single dish and interferometric data onto identical grids. 
The single-dish Fourier data is multiplied by the ratio of 
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the Point Spread Function (PSF) volumes of the inter- 
ferometric image to that of the single dish image. This 
scaled data is added to the Fourier interferometric data 
multiplied by a feathering function. The feathering func- 
tion used in AIPS++ is l-FT(PB dish ), where FT(PB dish ) 
denotes the Fourier transform of the primary beam of 
the single-dish. This combined data is then Fourier trans- 
formed back to the image domain. 

Although the angular resolution of the final image at 
327 MHz is 90" x 45" at a position angle of 164°, and 
the rms noise level is 10 mJy beam -1 , comparable to the 
earlier results of Dubncr ct al. (1991), the newly processed 
327 MHz image gives a more reliable representation of 
Puppis A at all spatial frequencies. 

2.2. New radio continuum observations at 1425 MHz 

We observed Puppis A at 1425 MHz using the VLA in 
the DnC and CnB configurations. The DnC data were 
obtained on 2003 January 21, 22, and 23, while the CnB 
observations were carried out on 2004 February 4, 5, 7, 
and 15. The flux density calibrators used were 3C 147 and 
3C 286 for the DnC observations and 3C 286 for the CnB 
observations. The phase and amplitude calibrator used for 
all the data observed was the source J0828— 375. Due to 
the large size of this SNR (about 55' in diameter) we used 
a mosaicking technique with 39 different pointings with a 
spacing of half of the primary beam width. Each sub-field 
was observed in two 50 MHz bands centered near 1465 
and 1385 MHz in continuum mode. The observations were 
made in cycles with 28 s in each pointing of the mosaic. 
This observing cycle was repeated 9 times in the DnC 
array and 12 times in the CnB configuration. 

The data were reduced using the AIPS++ software 
package. The data from each day were calibrated, imaged 
and self-calibrated separately to ensure that there were 
no day-to-day amplitude discrepancies. All the calibrated 
data sets were then concatenated into a single visibility 
datasct. We imaged and self-calibrated the combined vis- 
ibility data in an iterative fashion. The imaging was done 
using the standard mosaicking technique that jointly de- 
convolves all the fields, taking into account the primary 
beam effect at each pointing. The multiscale clean method 
was applied to deconvolve the image. We used uniformly 
weighted visibility data for the high resolution image. For 
the image that is used to compare with the 327 MHz data, 
the visibility data were tapered in such a fashion that the 
resulting restoring beam is identical to that of the 327 
MHz image. 

To recover the large scale structure missed by the in- 
terferometer, we used an image from the Parkes Southern 
Galactic Plane Survey (McClure-Griffiths et al. 2001). The 
single-dish and the interferometric mosaic images were 
combined using the same feathering technique described 
above. 

The uniform weighting image that was produced has 
a resolution of 34" x 16" at a position angle of —174°, 



and an rms noise of 0.5 mJy beam -1 , representing almost 
an order of magnitude improvement in sensitivity over the 
image of Dubncr ct al. (1991). 

3. Results 

3.1. The new VLA image of Puppis A 

Fig. 1 shows the new VLA 1425 MHz image of Puppis A. 
At this frequency the remnant has the appearance of an 
asymmetric shell with an average diameter of ~ 55' or 35 
pc at a distance of 2.2 kpc. The intensity of the brightest 
features along the periphery is spatially variable, ranging 
from ~ 8 mJy beam -1 in the N to ~ 60 mJy beam -1 to- 
wards the E. No sharp boundaries are observed, although 
the presence of a weak radio halo at a mean level of 3.5 
mJy beam -1 is apparent in the original images extending 
about 5' ahead of the bright rim. With the addition of 
short spacing contributions, all the flux density is repre- 
sented in the image in Fig. 1. This image shows diffuse 
interior emission at an average level of 8 mJy beam . 

These new sensitive radio observations show no evi- 
dence for any radio counterpart (3ct of 1.5 mJy beam -1 ) 
to the X-ray point source RX J0822— 4300 nor for a pul- 
sar wind nebula, consistent with previous negative results 
obtained by Gaensler et al. (2000). 

Enlargements showing the detailed 1425 MHz total in- 
tensity morphology in three selected regions of the Puppis 
A SNR are shown in Fig. 2. Due to the high dynamic 
range (w 1000:1) of the 1425 MHz image, several grayscale 
ranges are used to display the different relevant regions. 
Fig. 2a shows the E portion of the SNR shell, Fig. 2b 
displays the N part of Puppis A, while Fig. 2c shows 
the locations of the brightest features in the SW edge of 
the remnant. The new high resolution image reveals that 
the boundaries of Puppis A are quite structured, in con- 
trast with the smooth interior. The outer border consists 
of a series of arc-like features, resembling a "wave-like" 
structure. Such features are more evident in the N half 
of Puppis A (Fig. 2b) , although some are detected to the 
SW. We note that along the NW and SW borders, the 
"wave-like" features are observed tangential to the outer 
rim, but their orientation changes to approximately radial 
(i.e perpendicular to the local shock front) along the NE 
border. 

Petre et al. (1982) reported the detection of bright 
X-ray emission near R.A.= 08 h 24 m , Dccl.= -43° 00' 
(J2000), the BEK (see Fig. 2a). Intense radio emission 
is observed coincident with the X-ray emitting area (see 
Fig. 7 below). In addition, the new high resolution 1425 
MHz observations reveal the presence of a compact (size 
~ 45") feature of enhanced synchrotron emission (fea- 
ture A in Fig. 2a), exactly at the site where the shock 
front forms an angle, near R.A.= 08 h 24 m 18 s , Decl.= 
—43° 00' 00" (observed in Fig. 2a as a bright region ad- 
jacent to the BEK). To the S of this feature, the presence 
of a ~ 6' prominent, straight filament (labeled as B in 
Fig. 2a) is noticed near R.A.= 08 h 24 m 22 s , running N- 
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Fig. 1. The SNR Puppis A at 1425 MHz. This image was obtained from a combination of 39 different pointings observed with 
the VLA in the DnC and CnB configurations in 2004. Single dish observations from the Parkes Southern Galactic Plane Survey 
(McClure-Griffiths et al. 2001) have been added to the data. The angular resolution is 34" x 16" (beam shown at lower left), 
P.A.= — 174° with an rms of 0.5 mjy beam -1 . The brightness range covered by the grayscale is between 3 and 35 mjy beam -1 . 



S from Decl.= -43° 02' to Decl.= -43° 08'. To the E of 
this filament, a fainter semi-circular protrusion (extending 
about 6') is observed. The SW region of Puppis A (Fig. 2c) 
shows the only bright feature located on the edge of the 
SNR. Short arcs (size ~4'), similar to those observed to 
the NW, are also observed in this region, although in this 
case they appear displaced by ~ 8' from the limb to the 
interior. 

Fig. 3 shows the re-processed image at 327 MHz of 
Puppis A. These data were used with the 1425 MHz image 
for the spectral analysis discussed in Section 4. 

4. The spectral properties 

4.1. The integrated spectrum 

We have calculated integrated flux densities for Puppis 
A based on the new images, obtaining: S327 MHz =263 ± 
20 Jy and S 1425 MHz = 114 ± 8 Jy. The errors reflect the 



uncertainties both in the choice of background emission 
and integration boundaries, in addition to the inherent 
rms noise of the images. 



To construct accurate spectra of individual SNRs it is 
necessary that all flux density values be brought on to the 
same absolute flux density scale. Table 1 contains the flux 
densities between 19 and 8400 MHz. All the flux densities 
for frequencies above 400 MHz are on the Baars et al. 
(1977) scale. For those below that frequency no scaling 
was applied. Fig. 4 shows the resulting radio spectrum of 
the Puppis A SNR; the new determinations are plotted 
as filled circles. The flux density values are adequately 
fitted by a single power law of slope a = —0.52 ± 0.03, in 
good agreement with estimates presented by Milne & Hill 
(1969) (a ~-0.50) and Dubner et al. (1991) (a ~-0.53). 
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Fig. 2. Enlargements of the new VLA 1425 MHz image, (a) E region, containing the brightest features of the remnant. Feature 
A marks the position of a bright point-like emission region, while B indicates the location of the straight N-S filament. The 
grayscale ranges from 5 to 52 mjy beam -1 , (b) N part of Puppis A; the irregular outer border is prominent. Note the presence 
of "wave-like" features almost perpendicular to the shock front on the NE and parallel to it on the NW. The grayscale ranges 
from 6 to 16 mjy beam" 1 , (c) The SW region showing the presence of the only bright feature located at the outer boundary 
of Puppis A. The grayscale varies between 5 and 18 mjy beam -1 . 



4.2. The spatial spectral distribution 

It is important to discern whether the electron energy 
spectrum is a global property of this SNR, or if it is lo- 
cally affected by the neutron star and/or the shock-ISM 
interaction. Hence, we performed a careful study of the 
variations in the spectrum as a function of position within 
the remnant. 

The previous study of the spectral index distribution 
was carried out by Dubncr ct al. (1991) using VLA ob- 
servations at 327 and 1515 MHz. The distribution of the 
spectral indices was estimated by integrating flux densities 
within five rectangular boxes (of different sizes according 
to the region, varying between approximately 16 and 190 
arcmin 2 ), and comprising the brightest portions of Puppis 
A. An index as steep as a ~ —0.7 was suggested for the E 
border, while for the other four regions the spectral index 
was found to be about —0.5. Using more recent spectral 
analysis techniques and the new VLA 1425 MHz data, 
we are able to produce a detailed analysis of the spatial 



spectral transition from the bright features to the smooth 
weaker interior of this SNR. 



To accurately determine the spectral index distribu- 
tion from the comparison of two interferometric images 
at different frequencies, it is necessary to ensure that the 
corresponding images are matched in uv coverage. As men- 
tioned in Section 2.2, the data at 1425 MHz were processed 
by applying appropiatc weights to the visibilities to obtain 
an image with the same beam size and beam shape as the 
327 MHz image. Also, in order to avoid any positional 
offsets, the images were aligned and interpolated to iden- 
tical projections (field center, pixel separation, etc). In 
addition, we checked for possible flux density calibration 
bias. That is, an extrapolation back to zero flux density 
at one frequency should correspond to zero flux density 
at the other frequency, otherwise there is a zero-point off- 
set that must be corrected. We have checked the 327 and 
1425 MHz images, finding only an insignificant offset of 
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Fig. 3. Radio continuum image of SNR Puppis A at 327 MHz obtained with the VLA in the CnB configuration in f 988 and 
now re-processed. The interferometric data were combined with single dish data from MPIfR Bonn radiotelescope (at 408 MHz, 
Haslam et al. 1982). The final beam size is 90" x 45" (shown at lower left) at a position angle of 164°, and the rms noise level 
is 10 mjy beam" . The grayscale ranges from 35 to 350 mjy beam - . 



Table 1. Integrated flux densities for the SNR Puppis A 
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Fig. 4. Radio spectrum of Puppis A from the scaled flux den- 
sity values listed in Table 1. The data from this paper are 
shown by filled circles (327 and 1425 MHz). The linear fit to 
the values of flux density indicates a global spectral index of 
—0.52 ± 0.03. The dashed line shows the linear fit extrapolated 
to 19 and 86 MHz, at which frequencies no correction has been 
applied to bring the measured values to the Baars et al. (1977) 
scale. 



Frequency 


Scaled flux 


References 


(MHz) 


density (Jy) 




19 


800 ± 160 (a) 


Rishbeth (1958) 


86 


690 ± 100 (a) 


Mills et al. (1960) 


327 


263 ± 20 


This work (b) 


408 


235 ± 20 


Haslam et al. (1982) 


635 


180 ± 29 


Milne & Hill (1969) 


843 


144 ± 10 


Arendt et al. (1990) 


960 


130 ± 12 


Harris (1962) 


1410 


129 ± 20 


Milne & Hill (1969) 


1425 


114 ± 8 


This work (t) 


1440 


166 ± 17 


Mathewson et al. (1962) 


1515 


118 ± 10 


Dubner et al. (1991) 


2650 


92 ± 14 


Milne & Hill (1969) 


2700 


78 ± 12 


Milne (1971) 


4750 


59 ± 5 


Milne et al. (1993) 


5000 


67 ± 7 


Milne & Hill (1969) 


5000 


61 ± 7 


Milne (1971) 


8400 


38 ± 4 


Milne et al. (1993) 



No correction to Baars et al. (1977) scale was applied. 
' b ' Flux density scale from VLA Calibrator Manual, 
http:/www. aoc.nrao.edu/~gtaylor/calib. html. 



0.1 mjy beam 1 between frequency pairs. Therefore, no Based on these images, we carried out the spatial spec- 
further correction was applied. tral analysis. This was done using two different methods 
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in order to confirm that the observed spectral features are 
real, independent of any possible contamination from the 
databases used for the comparison. 

In Fig. 5 we show the spectral index distribution of 
Puppis A as obtained from the direct comparison of the 
matched images at 327 and 1425 MHz after clipping 
at the 5 a level for both frequencies. We find that the 
spectral index varies across the SNR approximately be- 
tween a ~ —0.8 and —0.2. In this grayscale representa- 
tion, darker means flatter spectrum. It is interesting to 
note that the feature with the flattest spectrum (seen as 
a black band running from NE to SW in the interior of 
Puppis A) includes the location of the neutron star RX 
J0822— 4300 and of the geometrical center of the explo- 
sion as derived from optical observations by Winkler ct al. 
(1988) (shown in Fig. 5 as a triangle and a plus sign re- 
spectively). These positions are known with an accuracy 
of tenths of an arcsec for RX J0822— 4300 and ~ 1 arcmin 
for the proposed center of the explosion. We can confi- 
dently assume that this spectral feature is real since both 
images have recovered all flux density. 

The present spectral map shows that along most of 
the periphery the spectrum is clearly steep (steeper than 
a=— 0.7). To the W, the ear-like feature shows a combina- 
tion of different structures with variable spectral indices. 

Curiously the spectral pattern observed near the 
bright eastern radio features (approximately near R.A.^ 

08 h 24 m , from Dccl. 42° 55' to 43° 10') does not 

correlate with the morphology observed in total power (see 
Fig. 1). Although it is somewhat flatter than the rest of 
the shell, it is apparent that the spectrum varies follow- 
ing a pattern of parallel fringes that alternate steeper and 
flatter spectral indices. 

In order to test the spectral index distribution de- 
termined above, we also constructed a gallery of tomo- 
graphic maps. Spatial tomography between two images is 
a technique developed by Katz-Stone & Rudnick (1997a), 
based on the use of a test spectral index a t - A tomo- 
graphic map is calculated by scaling the brightness in 
the higher frequency image (S1425MHZ, in our case) by 
the test spectral index, and subtracting this scaled im- 
age from the lower frequency image (8327 MHz)- Then, 
a tomographic gallery is obtained by calculating St = 
S327 MHz - (-Ms)"* s 1425 MHz, stepping through a range 
of spectral indices a t . Features which have a spectral in- 
dex identical to the test value will vanish in the tomo- 
graphic map. Spatial components that have different spec- 
tral indices will appear as positive or negative features 
depending upon whether the spectrum is steeper or flat- 
ter than the assumed test value. This method has been 
succcsfully applied to analyze spatial spectral variations 
in the SNRs of Tycho (Katz-Stone et al. 2000), Kepler 
(DeLaney et al. 2002), G292.0+1.8 (Gaensler & Wallace 
2003) and in the study of the radio galaxies 3C 67, 3C 190 
and 3C449 (Katz-Stone & Rudnick 1997a,b). The tech- 
nique is particularly useful to discntaglc components with 
different spectra that can overlap along the line of sight. 
Also, when the goal is to locate small-scale spectral vari- 




RIGHT ASCENSION (J2000) 

Fig. 5. Grayscale image of the spectral index distribution 
across the SNR Puppis A as derived between 327 and 1425 
MHz. The wedge displays the spectral index value (xlOOO). 
The 1425 MHz intensity contours plotted for reference are 36, 
60, 90 and 160 mjy beam - . A cutoff at 5 ct was applied at 
both frequencies. The plus sign indicates the geometrical cen- 
ter of the SN explosion, while the triangle marks the position 
of the neutron star RX J0822-4300. 



ations, this method can provide a more accurate picture 
than the simple comparison between images at different 
frequencies. 

In Fig. 6 we show the resulting tomographic images for 
four test values: a t =-0.2 (Fig. 6a), -0.4 (Fig. 6b), -0.6 
(Fig. 6c), and —0.8 (Fig. 6d). In these maps the spatial 
components appear as positive (light) or negative (dark) 
features for a spectrum steeper or flatter than the assumed 
test value respectively. The point sources overlapping the 
remnant are always seen as positive residuals in the tomo- 
graphic images (white spots) as they have spectral indices 
steeper than —0.8, thus supporting their extragalactic na- 
ture, as already noted by Milne et al. (1983) and Dubner 
et al. (1991). 

The most interesting tomographic images are those ob- 
tained from the comparison with at=— 0.4 (Fig. 6b) and 
with at=— 0.6 (Fig. 6c), where, as noticed before, a no- 
table striation is evident along the periphery. It is worth 
noting that this spectral pattern formed by short fringes 
with a alternatively steeper and flatter than the back- 
ground, mimics the "wave-like" morphology noticed in the 
total power image along the NE, NW and S borders. That 
is, the spectral components appear almost parallel to the 
outer boundary to the NW and SW, and transverse to it 
along the E side. 

The newly identified point-like feature to the E of 
Puppis A (feature A in Fig. 2a), has a spectrum steeper 
than —0.8 (i.e. it appears as a positive residual in all tomo- 
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Fig. 6. Series of tomographic images for the SNR Puppis A constructed from data at 327 and 1425 MHz, which were matched 
in uv coverage, (a) a t =— 0.2; (b) a t =— 0.4; (c) a t =— 0.6; and (d) a%=— 0.8. Bright areas have radio spectral indices steeper 
than the assumed at and dark areas imply flatter spectral indices. The same grayscale is used to display the four images. 

small scale structures that include filaments, short arcs 
and knots are, in general, poorly correlated with radio 
features. The brightest features are preferentially located 
to the E although the N-NW rim is also well delineated 
by bright X-ray filaments. Fig. 7 demonstrates that the 
most striking radio/X-ray correspondence occurs predom- 
inantly along the E limb of the remnant. Three features 
are well matched: the BEK, the curved filament near 
R.A.= 08 h 23 m 20 s , Decl.= -43° 02' (sec below) and the 

indentation around R.A.~ 08 h 24 m 30 s , Deck 43° 00'. 

In contrast, little radio/X-ray agreement is observed in the 
SE sector of Puppis A; the radio emission extends over 5' 
beyond the X-ray limb. Moreover, the protrusion ahead 
of the bright, straight N-S radio filament (feature B in 
Fig. 2a) lacks a counterpart in X-rays. A similar radio/X- 
ray comparison cannot be carried out to the SW of Puppis 
A due to a missing ROSAT pointing in this direction (see 
the ROSAT archive); thus the reconstruction of the im- 
age is likely incomplete in this region. Towards the NW 
border of Puppis A the X-ray emission accompanies the 
"wave-like" morphology revealed in the radio, with the 
same orientation parallel to the shock front. The dense 
X-ray knot located immediately behind the NW shock 
front (the "bright northern knot", BNK, R.A.~ 08 h 22 m , 
Decl.~ —42 ° 37' ) is prominent in both spectral domains. 
Petre et al. (1982) have suggested that this knot origi- 
nated as a result of the interaction with a dense neutral 



graphic images). We thus conclude that, in all probability, 
it is an extragalactic object. 

The central band seen in Fig. 5 disappears against the 
local background in the tomographic image traced with 
a t =— 0.2 (Fig. 6a), thus confirming its very flat spectrum. 

The grayscale corresponding to the at =—0.8 image 
(Fig. 6d) reveals that the spectral distribution in Puppis 
A has practically no components steeper than this value. 

5. Comparison between radio and X-ray emission 

The X-ray properties of Puppis A have been investigated 
using the Einstein, ASCA and ROSAT telescopes (Pavlov 
et al. 1999, and references therein). The highest resolu- 
tion X-ray image of the entire source was obtained with 
ROSAT (HRI) in the energy band from 0.1 to 2 kcV (Petre 
et al. 1996). 

In Fig. 7 we present a false color image showing the 
spatial comparison between the radio image at 1425 MHz 
(in red) and the soft X-ray emission (in green) for the 
entire remnant. For the first time a comparison with a 
spatial resolution < 30" is possible. Features where both 
spectral bands overlap are shown in yellow. 

The X-ray emission associated with Puppis A has been 
interpreted as resulting from several thermal components 
(Bcrthiaume et al. 1994). The N half of Puppis A is clearly 
the brighter and more complex in X-rays. A network of 
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Fig. 7. A high resolution X-ray/radio comparison of the Puppis A SNR. The green image corresponds to soft X-ray emission 
from ROSAT in the 0.1-2 keV range (Petre et al. 1996), while in red the 1425 MHz radio emission is shown. The yellow regions 
are areas where emission in both spectral ranges overlap. The X-ray image is smoothed to the resolution of the radio image, 
34" x 16", P.A --174°. 



cloudlet, similiar to the BEK. In the radio domain, the 
spectrum of this knot between 327 and 1425 MHz has a 
spectral index a ~ — 0.5 (see Section 4.2.2), similar to the 
global spectral index of this SNR. A careful spectral study 
of this clump in X-rays would help towards understanding 
its origin. 

Hwang et al. (2005) have used the Chandra X-ray tele- 
scope to resolve the fine structure around the BEK over a 
field of view of 8' x 8'. These authors describe the inner 
structure in this part of the remnant as consisting of: (1) a 
bright compact knot, (2) a curved vertical feature (named 
the "bar") about 1' W behind the shock, and (3) a smaller 
bright cloud (the "cap" ) separated from the "bar" by faint 
diffuse emission. Nomenclature for these X-ray features is 
contained in Fig. 8. 

To study the relation between the X-rays and the ra- 
dio components in the Puppis A SNR, in Fig. 8 we show 
a direct comparison of the new 1425 MHz VLA image 
with the Chandra X-ray image obtained from the combi- 
nation of the soft, medium, and hard X-rays (note that the 
Chandra image has an angular resolution about 5 times 



finer than the radio image). Fig. 8 reveals that the X-ray 
emitting plasma does not coincide with radio features at a 
detailed level. However, on a larger scale reasonable agree- 
ment is observed. The radio emission is compressed in 
the indented eastern X-ray boundary, confirming that the 
shock wave is wrapping around an obstacle in the neutral 
ISM as was earlier proposed by Hwang et al. (2005). A lo- 
cal radio maximum matches the BEK, as shown in Fig. 8. 
In this portion of the SNR, this feature is not the bright- 
est feature in radio as is the case for X-rays. The most 
intense radio emission is found to the E of the X-ray "bar" 
(around R.A.= 08 h 24 m 14 s , Decl.= -42° 59' 49"), where 
it correlates with a secondary X-ray peak (see Fig. 8). 
No conspicuous radio feature can be associated with the 
"cap" . The region of faint diffuse X-ray emission between 
the "bar" and the "cap" is a region of low surface bright- 
ness radio emission. 

In summary, based on the present comparison only 
moderate agreement between radio and X-ray emission 
has been observed at a detailed level. This conclusion is 
not surprising since the X-rays can be associated with 
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Fig. 8. (Left) ROSAT HRI mosaic of the Puppis A SNR. The box indicates the region of the E limb observed with Chandra 
(Hwang et al. 2005). The plus (+) symbol shows the position of the compact source RX J0822— 4300. (Right) The grayscale 
representation combines soft (~0.4-0.7 keV), medium (~0.7-1.2 keV), and hard band (>1.2 keV) images around the bright 
eastern knot, while the superposed contours trace the radio synchrotron emission at 1425 MHz at 18, 23, 29, 36, 43, 50, 56, 61, 
68, 74, 80, and 87 mjy beam -1 . 



thermal radiation from a hot plasma, while the radio con- 
tinuum arises via synchrotron emission from relativistic 
electrons. 

6. Is Puppis A a new case of a Galactic SNR 
shaped by energetic jets? 

As noticed before, the general shell-like morphology of 
Puppis A appears distorted to the E by the protrusion 

near R.A.~ 08 h 24 m 30 s , Dccl. 43° 06' and to the W 

by the "ear" around R.A.~ 08 h 20 m 00 s , Decl. 42° 54' 

(Fig. 1, Fig. 2a and Fig. 7). If we trace two intersecting 
lines oriented about a position angle of ~ 100° (from N 
through E), with a semi-opening angle of ~ 15°, crossing 
at the position of the neutron star candidate RX J0822- 
4300, the appearance of Puppis A strikingly resembles 
that of the SNR W50. The shell of W50 is distorted by 
the processing jets that originate in the central compact 
source SS433, and Fig. 9 illustrates the similarity between 
the two SNRs. In Fig. 9 (left) we have plotted a sketch of 
the hypothetical "cones" overlapping the new radio image 
of Puppis A, while on the right side an image of the SNR 
W50 at 1.4 GHz (from Dubner et al. 1998) is shown with 
the precession cones of ~ 20° semi-opening angle delimited 
by the intersecting white lines. In the case of W50, the 
elongated shape has been interpreted as being due to the 
action of the two oppositely directed relativistic jets in- 
jected by the X-ray binary SS433. The present comparison 
shows that not only the general appearance with symmet- 
rical lateral extension is similar for both SNRs, but that 
the analogy also extends to particular features. For exam- 



ple, the bright straight N-S filament in Puppis A (named 
B in Fig. 2a) looks very similar to the straight radio fila- 
ment that stands out near the E extremity of W50, and 
which is thought to be the radio "hot-spot" produced by 
the interaction of the relativistic jet from SS433 with the 
ISM (Safi-Harb & Ogelman 1997). Beyond these bright, 
straight filaments, diffuse radio emission is observed in 
both cases. 

The compact X-ray source RXJ 0822-4300 in Puppis 
A has been proposed to be the compact stellar remnant 
formed in the SN explosion (Pavlov et al. 1999). It is a 
radio-quiet neutron star belonging to the class of "Central 
Compact Objects" (CCOs) in SNRs whose nature remains 
unclear (Pavlov et al. 2002). Based on detailed HI ob- 
servations carried out in the direction of RXJ 0822-4300, 
Reynoso et al. (2003) suggested that this CCO could be 
ejecting two oppositely directed jets that are creating a 
bi-lobed cavity in the HI. This elongated hole is about 
15' long, oriented in the NE-SW direction, with a position 
angle of about 60°. Both the X-ray point source and the 
position of the geometrical center of the explosion (sug- 
gested to be at R.A.~ 08 h 22 m 27 s . 5, Decl. 42° 57' 29", 

Winkler et al. 1988) lie within the HI hole. Thus, the 
possible HI signature of jets is not exactly aligned with 
the direction of the hypothetical cones that connect the 
large scale E and W extensions. However, it cannot be dis- 
counted that the HI feature is related to the path of the 
neutron star from the center of explosion to its present 
position, while the large scale cones are likely to be asso- 
ciated with the persistent action of jets on the surrounding 
SNR shell. 
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Fig. 9. Comparison of the SNR Puppis A (left) with the SNR W50 that contains the central X-ray binary SS433 {right, from 
Dubner et al. 1998). In Puppis A, the white cross shows the location of the X-ray point source RX J0822-4300, a radio-quiet 
neutron star proposed to be the compact remnant of Puppis A explosion. The white intersecting lines delimit the "cones" 
proposed to be associated with Puppis A in analogy with the case of W50. In the image of the SNR W50, the white lines trace 
the volume defined by the oppositely directed, precessing, relativistic jets. 



In the X-ray domain, no conspicuous feature is appar- 
ent in the immediate vicinity of RXJ 0822-4300 (see Fig. 7 
and Fig. 8 left) . The first bright X-ray features within the 
hypothetical cones are the thin filament to the E near 

R.A.~ 08 h 24 m 15 s , Decl. 43° 05' and to the west near 

R.A.- 08 h 20 m 30 s , Decl. 43° 00' (see Fig. 8 left for a 

representation of the X-ray emission in equatorial coordi- 
nates). At the assumed distance of 2.2 kpc, these features 
occur at about ~ 32 pc E and ~ 22 pc W of the central 
point X-ray source. In the case of W50, the X-ray lobes are 
symctrically displaced E and W of SS433 starting at about 
15' from the center (approximately 17 pc at the distance 
of W50) and have a knotty structure with enhanced emis- 
sion at about 35' (i.e. about 40 pc) from SS433 (Safi-Harb 
& Ogelman 1997). The displacement of the X-ray lobes in 
W50 has been interpreted by Safi-Harb & Ogelman (1997) 
as due to the expansion of beams in a low density ISM. 

A careful study of the X-ray features in the interior 
of Puppis A is highly needed to understand the interac- 
tion of the central neutron star with the surrounding SNR 
and to answer the question as to whether Puppis A is an- 
other example of a Galactic SNR distorted by the action 
of symmetric energetic jets. 

7. Conclusions 

In the current work we have presented a new image of 
the SNR Puppis A at 1425 MHz. This image has revealed 
with superb detail the presence of short filaments along 
most of the distorted periphery of Puppis A. In particu- 
lar, in the northern half of Puppis A the radio emission 
near the limb resembles a "wave-like" structure, oriented 
parallel to the shock front along the NW but almost per- 
pendicular to it towards the brighter NE limb. Such a pat- 
tern is replicated in the spatial spectral distribution, where 
fringes with analogous orientation alternate between flat 
and steep spectra. We propose that the peripheral fila- 



ments are a manifestation of Raylcigh- Taylor (R-T) in- 
stabilities that arc distorting the interface between the 
ejecta and the ambient material. As R-T instabilities de- 
velop, they stretch and compress the magnetic field (Jun 
& Norman 1996). This produces magnetic field amplifica- 
tion in the ambient medium close to the interface, which 
in turn enhances the synchrotron emission. In our model, 
the R-T fingers observed to the NE are radial, while to 
the NW they are tangential suggesting that along this 
border the shock is probably in a more advanced evolu- 
tionary stage. This could be related to the possibility that 
to the NE the shock front is expanding in a denser ISM. 
To confirm this hypothesis, it is important to investigate 
the magnetic field orientation in the region. A detailed 
study of polarization vectors at high angular resolution is 
planned for the near future. 

Similar short, curved filaments are also observed in the 
young SNRs of Tycho and Cas A, both in the radio and 
in X-rays, although in these cases a weaker plateau and a 
thin brighter outer rim is observed ahead of the filaments 
(Reynoso et al. 1997; Hwang et al. 2004). In Puppis A, a 
faint plateau can be clearly seen around almost all the pe- 
riphery, but no sharp outer border, presumably associated 
with the foward shock, is detected in this intermediate-age 
SNR. 

The most conspicuous radio features are situated at 
the E border, at the site where the shock front has en- 
gulfed a dense ISM cloud. Interestingly enough, the radio 
counterpart of the BEK and the narrow straight filament 
located just south of it do not appear as spectrally dis- 
tinct features. If these synchrotron maxima are the result 
of the interaction with external clouds, then we are led to 
the conclusion that this condition is not reflected in the 
radio spectrum. The spectral map of Puppis A does not 
show diferences between the E and the W sides in spite 
of the gradient in density in the ISM (Petre et al. 1982). 
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Wc can therefore conclude that in this case the shocks are 
probably still non-radiative and the compression is not 
very high. This contention is supported by the lack of op- 
tical filaments in the boundaries and the generally poor 
radio/optical correlation (Dubner et al. 1991). 

Another spectral peculiarity in Puppis A is the fact 
that the spectrum steepens in general from the interior 
to the periphery. In the first-order Fermi, non-relativistic, 
test particle limit, stronger shocks yield natter indices up 
to the a=— 0.5 limit. Therefore, based on this acceleration 
model, higher emissivity regions should have flatter spec- 
trum, while steep spectra are expected in the more diffuse 
emission regions (see for example Anderson & Rudnick 
1993), exactly the opposite of what we find in Puppis A. 
This result is confirmed through the tomographic anal- 
ysis of the spectrum. It is interesting to note that the 
region of Puppis A with the flattest spectrum is a narrow 
band towards the interior, like a channel that connects 
the position of the neutron star RX J0822— 4300 with the 
presumed geometrical center of explosion. The continuous 
injection of relativistic particles by the neutron star could 
provide an explanation for this spectral characteristic. 

Based on a compilation of integrated flux densities over 
a wide spectral range, we have recalculated the global non- 
thermal radio spectrum of the SNR Puppis A. The fit 
produces a spectral index of a=— 0.52 ± 0.03 for the whole 
SNR. 

We have also used the new radio image to perform the 
first accurate radio/X-ray comparison for Puppis A. In X- 
rays, Puppis A is bright on the E half, and the emission 
seen by ROSAT between 0.1 and 2 keV is very structured. 
The new radio image, although much more sensitive than 
previous ones, does not show detailed correlation with the 
thermal plasma imaged in X-rays, except for the bright E 
features. Even there, the matching is not accurate when 
the radio emission is compared with the high resolution 
Chandra image. The only remarkable coincidence revealed 
by the new results occur along the N-NW boundaries, 
where the "wave-like" radio features near the shock front 
have their exact counterpart in the X-rays. A detailed X- 
ray study along the NW side, including the BNK, would 
be valuable. To the S-SW of Puppis A, we were unable to 
carry out a similar comparison because a ROSAT pointing 
is missing there. 

On the other hand, based on the new observations, we 
have confirmed that down to ~1.5 mJy beam -1 neither 
a radio counterpart to the neutron star RX J0822— 4300, 
nor a surrounding pulsar wind nebula, is observed. 

Finally, based on the radio images and their compar- 
ison with the X-ray emission, a striking similarity be- 
tween Puppis A and the SNR W50, which harbors the 
Galactic microquasar SS433 in its interior, is seen. This 
suggests that Puppis A is a candidate for being another 
SNR shaped by internal jets, probably at an earlier evolu- 
tionary stage. Detailed central X-ray images and spectral 
studies would be of great help to test this suggestion and 
to understand the nature of the central X-ray compact 
object. 
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